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ABSTRACT 

Host galaxies of highest redshift QSOs are of interest; they provide us with a valuable oppor- 
tunity to investigate physics relevant to the starburst-AGN connection at the earliest epoch of 
the Universe, with the most luminous black holes. 

Here we report an optical detection of an extended structure around a QSO at z=6.43 in 
deep z' and z r -band images of the Subaru/Suprime-Cam. Our target is CFHQSJ2329-0301 
(z=6.43), the highest redshift QSO currently known. We have carefully subtracted a PSF con- 
structed using nearby stars from the images. After the PSF (QSO) subtraction, a structure in 
the z'-band extends more than 4" on the sky (i? e =ll kpc), and thus, is well-resolved (16cr 
detection). The PSF-subtracted z r -band structure is in a similar shape to that in the z'-band, 
but less significant with a 3 a detection. In the z'-band, a radial profile of the QSO+host shows 
a clear excess over that of the averaged PSF in 0.8-3" radius. 

Since the z'-band includes a Lya emission at z=6.43, we suggest the z' flux is a mixture 
of the host (continuum light) and its Lya emission, whereas the z, -band flux is from the 
host. Through a SED modeling, we estimate 40% of the PSF-subtracted z'-band light is from 
the host (continuum) and 60% is from Lya emission. The absolute magnitude of the host is 
Mi4 50 =-23.9 (c.f. Mi45o=-26.4 for the QSO). A lower limit of the SFR(Lya) is 1.6 M yr" 1 
with stellar mass ranging 6.2 x 10 s to l.lxlO lo M when 100 Myrs of age is assumed. The 
detection shows that a luminous QSO is already harbored by a large, star-forming galaxy in 
the early Universe only after ^840 Myr after the big bang. The host may be a forming giant 
galaxy, co-evolving with a super massive black hole. 

Key words: quasars:individual, cosmology : early universe, black hole physics, galaxies: high- 
redshift 



1 INTRODUCTION 

There have been accumulating evidence that galaxies and ac- 
tive galactic nuclei (AGN) co-exist. Evidence o f starburst fea- 
tures have been detected in n umerous AGNs ( Maioli no et al.l 
1 19971 : ICid Fernandes et alj l2005t) . Half of the ultraluminous in- 
frared galaxies contain simultaneously an AGN and starburst ac- 



trared galaxies contain simultaneously an ALrJN and starburst ac- 
tivity JGenzel et alj 1 19981 ; iGotol l2005h. P ost- starburst signatures 



have been found in AGNs dBrotherton et al.fl999l : lDewangan et al.1 
l200d : lGotd2006h . The discovery of tight correlation between black 
hole mass and bulge velocity dispersion provides another evidence 
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that the formati on of bulge and the central black hole may be 



closely linked dMagorrian et al J 1 19981; iFerrarese &~M erriti 



Gebhardt et al. 2000; Granato et al. 2001; Heckman et al 



2000 



2004 



Onken et alJl2004l) .The correlation implied the co-evolution of a 
super-massive black hole and its host galaxy. To understand galaxy 
formation theory and the AGN mechanisms, it is important to char- 
acterize a possible relationship between the galaxy formation and 
AGNs. 

At high luminosity end of the AGN, QSOs are 
not exception; At low redshift, QSOs are predominantly 
hosted by luminous, massive, bulge dominated galaxies 
dHamilton. Casertano. & Turnshekl |2002l: bunion et al.1 120031 : 



IPagani. Falomo. & Treves! 20031 : iFlovd et alj " 2004 ). irrespective 
of radio power. At high redshift, evolution of t he QS O host 
has been reported; ISchrammTw isotzki. & Jahnkej {2008 ) found 
hosts with stellar populations indicating recent massive star 
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Table 2. Basic properties of the constructed PSFs. The pixel scale is 
(0.2"/pix). The FWHM is obtained through a Gaussian fit. 



Figure 1. Composite pseudo-color image of the QSO (CFHQSJ2329- 
0301). The RGB colors are assigned to z',z r and i'-bands, respectively. 
The figures are north up, east left. 



formation at z ~ 3 . Aretxaga. Boyle. & Terlevichl dl995h : 
I Aretxaga. Terlevich. & Bovld i ll 9981) found SFR of >100- 
200M©yr~ 1 for z~2 QSO hosts. Ivillforth. Heidt. & Nilssonl 
d2008h als o found moderate rate of SFR for QSO hosts at z=0.87 
and 2.85. iHutchings etaD d2002L l200rj) found that hi gh redshift 
QSO host galaxies are irregular and highly disturbed. Hutchings 
d2005h reported a significant young stellar population with a 
possible spiral structure in two hosts at z=5. However also note 
that some luminous QSOs are reported to reside in luminou s 
massive early-type gal axies at 1 < z < 2 jKotilainen etai]|2007h . 
and at 2 < z < 3 ( Falomo et ai] |2008h . The host galaxies of 
radio-quiet QSOs are on average less luminous than the ones o f 
radio-loud QSOs (e.g jKukula et al.ll200ll : [Kotilainen et alj|2007l) . 
The apparent co-evolution of QSO and host galaxies leads to the 
question if QSO activity and star-formation are caused by the same 
physical mechanism, or one triggers the other. 

Although a proper investigation of high-z QSO hosts would 
shed light on the QSO-galaxy co-evolution, QSO host galaxies at 
highest redshifts of z > 6 have not yet been studied extensively, 
since due to high QSO-to-host light ratios, extremely deep, high 
resolution images are required to assure a sufficien t ly h igh signal- 
to-noise ratio. A pioneering attempt by Hutchings] d2005h did not 
resolve host galaxies of QSOs at z=6.23 and z=6.35. 

In this paper, we investigate a host galaxy of a QSO at z=6.43 
using deep optical images taken with the Suprime-Cam/Subaru 
equipped with new red-sensitive CCDs. This is the highest red- 
shift QSO as of the date of writing, and thus, provides us with 
an unprecedented opportunity to shed light on the co-evolution of 
galaxy and QSO at the massive end of super massive black holes 
with ~ 1O 9 M0, at the earliest epoch when the age of the Universe 
was only ~840 Myrs old. 

Unless otherwise stated, we adopt the WMA P cosmology: 
(h, fi m , Q L ) = (0.7, 0.3, 0.7) jKomatsu et alj2008h . 



Filter 


FWHM (pix) 


ellipticity 


i' (A e =768nm) 


2.45 (=0.49") 


0.04 


z' (A e =911nm) 


2.25 (=0.45") 


0.06 


z r (A e =988nm) 


2.69 (=0.54") 


0.07 



2 OBSERVATION 

Our target is CFHQSJ2329-0301, which is a QSO at z=6.43 
JWillott et al]|2007l) . Characteristics of the target is summarized in 
TableQ] Assuming Eddingt on luminosity, it s black hole mass is es- 
timated to be l.OxlO 9 M dFan et alj|200ll) . 

Observations were carried out on August 1, 26, and 27, 2008 
with the Suprime-Cam dMivazaki et al.ll2002f) . on which new red 
sensitive CCDs have been just installed (Utsumi et al. in prep.). 
The sensitivity of the new CCDs is 1.3 times better in z'-band, and 
1.9 times better in z r -band (A e // eC i it)e =988 nm). The data were 
originally taken to investigate the QSO environment, which will 
be presented elsewhere (Utsumi et al. in prep.). Since the observa- 
tion was partly for an engineering purpose to test the new CCDs, 
we used various exposure times with 360s x 10 in i', lOOsec x 7, 
300sec x 2, 400sec x 3, 500sec x 6, 700sec x 2 in z and 500sec 
x 19 in z r . The dithering pattarns were mostly circular with sepa- 
rations of 0.5-2 arcmins. The total exposure time on the target was 
3600 sec in i' , 6900 sec in z', 9500 sec in z,-bands. We used a 
nearby field observed both by the SDSS and ourselves for a pho- 
tometric calibration. A slight difference bet ween the SDSS and the 
Subaru filters are compensated using the iGunn&StrvkeiU 19831) 
stellar spectral library. The data are reduced using the modified 
version of the pipeline. The details of the reduction is described in 
Utsumi et al. (in prep.). The 2a limiting magnitudes are i'=26.73, 
z'=25.79, and z r =25.09 ABmag within 2" aperture. Here, la sky 
magnitudes are computed by randomly placing a 2"aperture in th e 
blank (sky) part of the image jYagi et alj2002l:IOuchi et alj2004l) . 

The pixel scale of the Suprime-Cam is 0.2'7pixel. The see- 
ing size of the images are ~0.5 arcsec (see Tabl^Ul. An i , z and z r 
composite image of the QSO is shown in Fig[T] where the RGB col- 
ors are assigned to z' , z r and i'-bands, respectively. Without even 
subtracting the central PSF, it is immediately noted that the QSO 
has an extended structure in the z'-band (red in the composite). 



3 RESULTS 

3.1 PSF subtraction 

In this section, we subtract a PSF from the QSO image to investi- 
gate a possibly extended structure around the QSO. In this process, 
it is critical to carefully construct the PSF model using the observed 
data since at such a high redshift the PSF has to be known very pre- 
cise to be able to see a marginally extended host galaxy. Using ob- 
served stars instead of analytical functions is especially important 
for z r -band where 0.5% of the flux is known to bleed in horizon- 
tal spikes due to the filter characteristics. In each band, we have 
carefully picked ~20 nearby stars within 200" in the mosaic image 
to avoid a possible PSF variation across the image. These stars are 
picked to be relatively isolated, and ambiguous stars according to 
their radial profile. Stars blended with a nearby object are excluded 
from the selection. We have specifically picked stars as bright as, 
but not brighter than the QSO itself to avoid a possible non-linearity 
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Figure 2. The upper, middle, lower panels are for %' , z' and 2 r -band, respectively.In each line, the left panels are reduced images. The middle panels are PSFs 
constructed using nearby stars. The right panels show residuals from the PSF subtraction. All figures are north up, east left. 



Table 1. Target information adopted from Willo tt et all fe007t) 



Object 


RA 


DEC z 


i' AB z' AB J M1450 


CFHQS J232908-03O158 


23:29:08.28 


-03:01:58.8 6.43 


>25.08 21.76±0.05 21.56±0.25 -26.4 



effect at the central pixels of bright stars. We have checked a choice 
of a different set of PSF stars does not affect results. Then we used 
DAOPHOT in iraf to construct an average PSF. The constructed 
PSF models in each band are shown in the middle panels of Fig(2] 
Basic properties of the PSF in each band are presented in Table [2] 



structed PSFs from the QSO images shown in the left panels of 
Fig|2] At first, gal fit was set to fit only the PSF and any remain- 
ing sky. The residuals are shown in the right panels of Fig[2] We 
discuss results from each band in the following subsections. 



Then we use galf it dPeng et al.|[2002h to subtract the con- 
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Table 3. Magnitudes and results of the fit. 



Object 


i' 
AB 


z ' 
AB 


ZrAB 


QSO+host 


25.54±0.02 


21.165 ±0.003 


21.683±0.007 


Host 


>25.34 (lcr limit) 


23.5±0.3 (16<t) 


24.3±0.2 (3cr) 


lcr sky 


25.44 (26 pix diameter) 


24.90 (26 pix diameter) 


25.46 (18 pix diameter) 



3.1.1 i' -band 

In the i! band (the top panels), no flux is left on the residual image, 
i.e., the QSO is consistent with a PSF in the i' band. This is not 
surprising since bluewards of the Lyman break, little light from the 
host galaxy can escape from the heavy absorption by the neutral 
hydrogen. 



3.1.2 z' -band 

In the z'-band shown in the middle panels of the Fig|2] there clearly 
remains an extended structure, subtending 4" after the PSF subtrac- 
tion. There is a little over- subtraction at the very centre, making the 
structure looks like a dough-nut shape perhaps because the PSF fit 
is affected by the extended structure. We have applied the same 
PSF subtraction to nearby stars for consistency check in Appendix 
lAl None of the PSF fits to stars yielded such an extended struc- 
ture in Appendix [A] We performed an aperture photometry using a 
13 pix (=2.6") of radius on the structure. The measured magnitude 
is z'=23.5±0.3. We have computed a l-a of the sky by randomly 
placing the same size aperture on the blank regions of the image. 
The 1 a sky noise is estimated to be 24.90 in z' band. Compared 
to the sky noise, the detection is 16u, i.e., the structure is well- 
resolved, and detected. 

In addition to the PSF subtraction, in the left panel of Fig[3] 
we show z'-band radial profiles of QSO+host in the blue solid line, 
and of the constructed PSF in the red dashed line. The profiles are 
normalized to be unity at a maximum. The solid, dashed, and short- 
dashed lines connect median in each sample. The error bars denote 
the errors in determining the median at each radii. At above the 
radius of 4 pixel (>0.8"), the figure shows that the QSO+host is 
clearly more extended than the PSF, showing the extended structure 
is well resolved. 

We note that there may be a point-source like structure to 
the right of the residual image in the z'-band in the middle-right 
panel of Fig[2] However, this is not likely a contamination from a 
foreground star since the possible point-source is not detected in 
i'-band, and thus quite red. Instead, it is more likely to be a star- 
forming region/substructure associated with the QSO. 



3.1.3 z r -band 

In the z r band image in the bottom left panel of Fig[2] there may 
be a possible elongation from north-east to south-west, in a similar 
orientation to the extended structure in the z'-band. Although the 
residual image on the bottom-right panel is much less obvious than 
in the z'-band, there is a hint of remaining flux. We performed an 
aperture photometry with an optimal 18 pixel diameter, obtaining 
24.3±0.2 mag of flux. The la sky noise measured by randomly 
placing the same size apertures on the blank regions of the image is 
25.46 mag, i.e., the extended structure is at a 3 a level. The errors 
on the PSF subtraction is 0.2 mag, and thus smaller. We also run the 




Figure 4. Both panels show residuals from the PSF subtraction in the z r - 
band. The right panel is box-car smoothed with lOpix. The figures are north 
up, east left. 

source extractor on the PSF subtracted image, detecting the host at 
a 2.6 a level over the local background. 

The radial profiles in the z r are shown in the right panel 
of Fig(3] Although the extended structure is less obvious, the 
QSO+host has a la excess over the PSF in three consecutive radii 
at 6-9 pixel. 

The extended structure becomes more visible in a 10-pix box- 
car smoothed image in Fig(4] For a sanity check, we fit the PSF to 
nearby stars in Appendix lAl where we do not see any clear excess 
in the residuals even after box-car smoothed. 

Because the Lya emission at z=6.43 is redshifted into the z'- 
band wavelength range, the extended structure in the z'-band could 
be an extended nebula of Lya emission (Lya blob) extensively 
searched at z~3 jMatsuda et alj|2004l ; lsaito et alj|2008h . However, 
in this QSO image, the extended structure is also found in the z r - 
band where the Lya emission does not contaminate the flux. There- 
fore, it is more likely that at least some z'-flux stems from the host 
galaxy (continuum), and the z'-band flux is a mixture of host and 
its Lya emission. We estimate the relative strength of the flux of 
the host (continuum) and Lya emission in detail in Section [X2l 

3.2 SED, stellar mass, SFR 

In Fig[5] we show the SED of the QSO (the black filled circles) 
and the host (the red rectangles). Since the host was not detected in 
the i', 1 a upper limits are shown at the effective wavelength. We 
computed the 1 a of the sky noise by randomly placing an aperture 
of the same si ze on the blank re g ion of the images, following the 
prescription bv lYagi et al.1 d2002l) : IOuchi et all fetXMl) . 

Overlaid are constant s tar-formation and delta burst mod- 
els ( iBruzual & Charlotl r2003h with the ages of lOOMyrs scaled 
to the z r -band flux. Assumed are the Salpeter initial mass func- 
tion, solar metallicity and no dust extinction. Flux bluewards 
the Lya emission is a t tenuat ed using the prescription given by 
iFan, Carilli. & Keatind fe006h . These models are often used to 
compute colors of Lyman break galaxies at z ~6. 
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Figure 3. Radial profiles of QSO+host (blue solid line), the constructed PSF (red dashed line), and the PSF+Sersic model (green short-dashed line) in the 
z'-band (left). The right panel is for the z r -band. Profiles are normalized at a maximum value. The pixel scale is 0.2"/pix. 



Comparing the colors of the host to the galaxy models, one 
immediately notices that the z'-band flux is not consistent with 
the flat continuum of the SED models, regardless of the star for- 
mation history. Perhaps this is because t he host galaxy has s ome 
Lya emission in the z'-band. In the past, iDavies et alj d2008l) de- 
tected extended line emission in high redshift galaxies on scales 
of ~4kpc. Such extended Lya emission neb ulae called Lya blobs 
have been exten sively searched at z ~3 dMatsuda et alj [20041 ; 
ISaito et alj[2008h . A QSO BR 1202-0725 has a Lya emitting com- 
panion, possibly ionized by the QSO |Hu. McMahon. & Egamil 
I l996t lOhvama. Taniguchi. & Shiovdl2004l). A typical FWHMo f 
Lya blobs is 300 km s" 1 dMatsuda et al.ll2004l; ISaito et al.ll2008h . 
If the Lya emission around the QSO is extended, it may be the- 
oretically predi cted infallin g gas to form a galaxy, illuminated by 
a central QSO dReesl ll988). If so, quantifying this will provide us 
with an important constraint on the formation of galaxies such as 
gas fraction, AGN geometry, and covering factors in the early Uni- 



Both of the SED models show an almost flat continuum slope 
at this small wavelength range between z and z r -bands. Thus scal- 
ing from the z r -bands flux and the SED models, in the z'-band, 
~2.5xl0 -19 erg cm" 2 s" 1 A" 1 flux is estimated to be the light 
from the continuum of the host galaxy. Subtracting this from the 
observed flux in the z'-band (6.0xl0" 19 erg cm" 2 s" 1 A" 1 ), 
~3.5xl0 -19 erg cm -2 s _1 A -1 of flux is estimated to originate 
from Lya emission of the host galaxy, i.e., in the z'-band, the con- 
tinuum to Lya flux ratio is 4:6. 

If such flux came through the Lya emission of 300 km s - , 
then the Lya luminosity would be 1.6xl0 42 er g s" 1 , which are 
comparable luminosity t o Lya emitters at z ~6 ( Kashikaw a et alj 
2004 lOuchi et alj|2008t) and Lya blobs at z=3.1(Matsu da et al.1 
20061) . 



Based on t he Lya luminos ity, we estimate the SFR(Lya) is 
> 1 .6 Mq yr - 1 dKennicuttl 1 9 98 ). However, the SFR estimated here 
are lower limits because significant fraction of bluer side of the Lya 
emission could be absorbed by the neutral hydrogen. 

We can also estimate star f ormation rate of t he host galaxy 
based on the UV luminosity using lKennicuttl dl998l) . Since the UV 



continuum of the host is relatively bright Mi45o=-23.9, we ob- 
tain SFR~50 Mq yr -1 . It has been suggested that SFR derived 
from UV continuu m is often larger than that from Lya, especially 
for L ya emitters dHu et alj|2004l ; lAiiki etai1l2003ljKodaira et all 
2003), as is the case here. Even though, this is quite a large number 
compared to that obtained through Lya (> 1.6 Mq yr -1 ), and thus, 
suggests the nature of the host galaxy; its UV continuum is brighter 
than typical Lya emitters. The host may be a massive galaxy that 
is decreasing its star formation activity. 

In the literature. Ijahnke et al.1 d2004h reported 2-30 Mq yr -1 
for z=l. 8-2.6 QSO hosts. IVillforth, Heidt. & Nilssonl d2008h mea- 
sured 0.03-33 Mq yr" 1 for QSO hosts at z=0.87 and z=2.75. The 
SFR estimated by UV light, however, is also subject to reddening 
uncertainties since a possible presence of dust has not been taken 
into account. 

With the lack of infrared photometry, it is difficult to constrain 
age/stellar mass through the SED. However, in the literature, a very 
young age of 5-100 Myrs for high re d shift galaxies a t 5< z <7 
has been reported dEgamietal.ll2005l ; lLai et alj[2007h . If we as- 
sume an age of 100 Myrs and scale models to z, -band flux, we 
obtain a stellar mass of l.lxlO lo M0 for the delta burst model, 
and 6.2 x 10 s Mq for the constant SFR model. Considering the age 
of the Universe is 840 Myrs, the host may be a forming massive 
galaxy, possibly co-evolving with the super massive black hole. 



3.3 Sersic fit 

Since we have a clear detection of extended structure in the z'-band, 
we attempted to fit the z'-band image with a PSF+sersic+sky model 
using the galfit. The z r -band image does not have enough 
signal-to-noise ratio to include Sersic profile in the fit. Results are 
shown in Fig[6] The middle panel of Fig|6] shows the PSF+sersic 
(PSF convolved) model, which is subtracted from the original im- 
age shown in the left panel. Some residual can be seen in the right 
panel, possibly i mplying an irregular m orphology of the host, being 
consistent with Hutchings et al. (2006) who pointed out that high 
redshift QSO host galaxies are found to be irregular and highly dis- 
turbed. Measured sersic parameters are _R e =2.0"±0. 1, axis ratio of 



Goto et al. 





I 



Figure 6. (left) The original z'-band image of the QSO. (middle) the best fit PSF+sersic model, (right) residual from the model subtraction. 
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who observed in J and Ks, found a mean of 11.1 ± 5.7 
kpc for their z w 1 sample. For their z ~ 1.5-2 sample they 
found a value of 9.78 ± 6.58 kpc for objects where it could be 
determined and upper limits (<10 kpc) for the remaining sources. 
iKotilainen et al.l (2007) found valu es of 6.7±1.7 kpc fo r their z fa 
1-2 sample (observed in K and H). lFalomo et al. (2008) found val- 
ues between 2.6 kpc and 11.3 kpc for a sample of three high redshif t 
QSOs (z « 2-3) observed in Ks. lVillforth. Heidt. & Nilssoj d2008t) 
found 1.5-3.9 kpc at z=0.9 and 1.5-3.9 kpc at z=2.8 in the FORS 
deep field. 

However, direct interpretation of our results is difficult; in the 
last subsection we estimated that the z'-band flux is a mixture of 
40% continuum light and 60%Lya emission. Therefore, this ex- 
tension of 11 kpc scale can be either extended Lya emission or 
the extended stellar population of the host. Clear morphological 
separation of these two have to await high-resolution near-infrared 
imaging, or specially-resolved IFU spectroscopy in the red end of 
the optical CCD. A spectro-polarization is also important to con- 
strain contribution from dust scattered light (Young et al. 2009, 
arXiv:0905.4l02l. 



Figure 5. SEDs of QSO and its host galaxy. Overplotted are SED models 
of constant SFR and delta starburst with 100 Myrs of age. The host is not 
detected in i'-band, where la upper limit is shown. 

b/a=0.74±0.03, and n = 0.019, suggesting that the host possesses 
a flatter profile than the PSF, and even than de Vaucouleurs. In the 
left panel of Fig(3] we show the PSF+Sersic model in the green 
short-dashed line. The PSF+Sersic model in the short-dashed line 
is a good fit within 8 pixels of radius, although the model underes- 
timates the profile significantly at >8 pixels of radius. 

The extended structure has an effective radius of i? e =2.0"± 
according to the sersic fit. Even in the image, one can easily rec- 
ognize the light distribution beyond 13 pixel (=2.6")- At the dis- 
tance of z=6.43, 2.0" corresponds to 1 1 kpc. As a reference, in the 
GOODS fields, a mean size of z ~ 5 galaxies wit h 1-51/* is 0.3" 
(~2 kpc) in restframe 1500A dFerguson et al.llioolh . In the UDF, a 
mean size of z ~ 6 galaxies with 0.1 -1.0L* is 0.15" (~1 kpc) in 
restframe 1600A( lBouwens et alj|2004l) . 

On the size of QSO hosts in the literature, iKukula et al.l 



4 DISCUSSION 

In Section|3] we detected an extended light in the z and z r -band, 
and suggested that this was the QSO host galaxy at z=6.43. How- 
ever, it is always a possibility that the extended light could be from 
a foreground lensing galaxy, which also explains the bright mag- 
nitude of the QSO by lensing magnification. However, we believe 
it is not very likely in this case; using the la limit of i' >25.34, 
the i — z 1 color of the host is > 2.1, which is very red and is 
hard to be explained unless the host has a break between i' and z . 
One possibility is a Balmer break at z=1.25. However, the color of 
the B aimer break is bluer, at approxima tely 0.64< u — g <1.99 
dFukugita. Shimasaku. & Ichik awalll995h . The size, absolute mag- 
nitude, nuclear-to-host flux ratio all do not seem to be unreasonable 
for the structure to be at z=6.43. In addition, no significant flux is 
detected at 6500A< A <8900A in the GM OS spectrum with one 
hour of exposure time jWillottetalJl2007l) . Therefore, it is more 
plausible to consider the extended emission is from the host galaxy 
of the QSO at z=6.43. 

We have estimated the host luminosity to be Mi4so=-23.9, 
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i.e., the QSO to host flux ratio in restframe UV is 8.8. This 
is a h igher ratio compared to the previous work. IHutchingsl 
d2005h found nuclear-to-extended ratio of 2.3 in rest 2100A. 
iGuvon. Sanders. & Stocktor] d2006l) report a flux ratio of 2.5 
(note their pass-bands are in NIR). If t he mass of the black 
hole is proportional to that of the host dMagorrian et al.l [l998l ; 
llvanov & Alonso-Herrero|[2003h . and if the QSO emits at a fixed 
fraction of the Eddington luminosity, then one would expect a cor- 
relation between the luminosities of the host and QSO. However, a 
number of reason can explain the deviation; the host may be more 
actively star-forming at earlier epoch; the Eddington ratio may be 
differe nt at z=6.43; the amoun t of nuclear obscuration may be dif- 
ferent. Hutchin gs et al. (2006) reported a higher QSO-to-host flux 
ratio for more luminous QSOs among their sample. Since this is 
only one example based on the restframe UV photometry, it is im- 
portant to obtain infrared photometry to accurately measure the 
stellar mass/ luminosity ratio of more hosts. High resolution spec- 
troscopy is also important to examine M — a relation at z>6. 

Comparing the m orphology of the host to those in the litera- 
ture, IHutchingsl J2003h resolved the host galaxies of 4 high redshift 
QSOs (z ~ 4.7) and found heavily disturbed galaxies, partially 
showi ng signs of ongoing merging. IGuvon. Sanders. & Stocktonl 
(2006) reported that 30% of PG QSO host galaxies showed ob- 
vious sign of disturbances, with strongly disturbed hosts favoring 
more luminous QSOs. In this work, there were some residuals after 
the host was fit with a sersic profile in the right panel of Fig|6] 
suggesting a possible non-smooth morphology. Further morpho- 
logical investigation, however, is difficult due to the limited depth 
and spatial resolution of our ground-based data. It is necessary to 
take high-resolution near-infrared AO or HST data to investigate 
the morphology of the host galaxy. 



5 SUMMARY 

We have detected an extended structure around a QSO at z=6.43 in 
the z and z r -band images of the Subaru/Suprime-Cam. The QSO 
is the highest redshift QSO currently known. After a careful PSF 
subtraction, the structure is detected at a 16 a level in z'-band and 
at a 3a level in the z r -band. In the z'-band, the radial profile of the 
QSO+host shows clear excess at > 4 pix of radius over the care- 
fully constructed PSF. There is a marginal excess in the z r -band 
profile at 6-10 pix of radius as well. Through the SED modeling, 
we estimate 40% of the PSF-subtracted z'-band light is from the 
host (continuum) and 60% is from Lya emission of the galaxy. 

The host is a bright, large galaxy with Mi45o=-23.9. Its ex- 
tended structure has an effective radius of >11 kpc. The lower 
limit on the SFR(Lya) is >1.6 Mq yr~\ the stellar mass of the 
host ranges from 6.2 x 10 s to 1.1 X 1O 1O M0 depending on the star- 
formation history when 100 Myrs of age is assumed. Therefore, 
the host is a large, star-forming galaxy, with a moderate amount of 
stellar mass. 

These results present an important example that indeed a super 
massive black hole resides in a large galaxy even in the early epoch 
when the Universe is ~840 Myr old. Quantifying the host mass 
and fraction of ionized gas in more detail will provide us with an 
important constraint on how galaxies and AGN form in the early 
Universe, such as gas fraction, AGN geometry, and covering fac- 
tors. 
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Figure Al. Reliability of the PSF subtraction is tested by fitting the con- 
structed PSF from 5 nearby stars, in i' and z'-bands. In each band, the left 
panels show actual stellar images. The right panels show residuals from the 
PSF-subtraction. 



APPENDIX A: PSF SUBTRACTION TEST 

Fitting the constructed PSF (Sec. 13. il l to the observed stars around 
QSO provides a sanity check of the PSF subtraction. In this exer- 
cise, we first select stars within 200" around the QSO as objects 
with source extractor's stellarity of > 0.97. We have specifically 
picked stars as bright as, but not brighter than the QSO itself to 
avoid a possible non-linearity effect at the central pixels of bright 
stars. These are the stars used to construct a composite PSF in 
Section [3] In Fig lAll we show 5 images of stars (left panels) and 
PSF-subtracted images (right panels) in i'(left) and jz'(right)-bands. 
Fig |A2l shows the same PSF subtraction for jz r -band, but the right- 
most panels show 10-pixel box-car smoothed images of residuals 
(the middle panels) in the z,--band. 

None of the residuals in Fig lAll shows an extended structure 
such as seen in the z' in the middle right panel of Fig[2] In the 
2 r -band, no obvious residuals are seen even after the images are 
box-car smoothed (c.f. Fig|4j. These sanity checks suggest that the 
extended structures found in Figs.|2]and|4]are not artifacts from the 
PSF subtraction. 



